Hybrid noiseless subsystems for quantum communication over optical fibers 
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We derive the general structure of noiseless subsystems for optical radiation contained in a se- 
quence of pulses undergoing collective depolarization in an optical fiber. This result is used to 
identify optimal ways to implement quantum communication over a collectively depolarizing chan- 
nel, which in general combine various degrees of freedom, such as polarization and phase, into joint 
hybrid schemes for protecting quantum coherence. 
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The ability to maintain and transmit coherent quan- 
tum superpositions is a prerequisite for the realization of 
protocols in quantum information processing. In commu- 
nication tasks, such as quantum key distribution |j| and 
teleportation |2| , a natural choice for the physical carrier 
of quantum information is optical radiation, often con- 
fined to optical fibers in order to minimize uncontrolled 
interactions with the environment. Quantum states en- 
coded in low-intensity, narrowband optical fields experi- 
ence two dominant decoherence mechanisms when trans- 
mitted through optical fibers [lj. The first mechanism 
is a linear loss of the light amplitude through residual 
absorption, scattering and coupling out to non-guided 
modes. The second one is a random transformation of 
the polarization state of the transmitted light due to 
changes in fiber birefringence caused by environmental 
conditions such as temperature fluctuations and aeolian 
vibrations. In this paper we address the latter mechanism 
and analyze ways to protect quantum coherence against 
depolarization by using general, multiphoton states of 
optical radiation. We assume here that linear losses are 
polarization-independent and that they can be dealt with 
in a standard way by postselecting transmissions when no 
photons have been lost. 

What makes fiber depolarization manageable in quan- 
tum communication tasks is that its temporal variations 
are relatively slow. In practice, the fiber birefringence re- 
mains virtually constant for long sequences of pulses that 
can be temporally resolved by the receiver, or even for 
roundtrip transmissions from the receiver to the sender 
and back. This enables a number of strategies for imple- 
menting communication protocols over such collectively 
depolarizing channels. The first one is to use a degree 
of freedom that is not affected by depolarization at all, 
such as the relative phase between consecutive pulses. 
This idea is best illustrated by time-bin encoding used 
successfully in a variety of experiments 0- A practi- 
cal issue in phase-encoding is the necessity to maintain 
a suitable phase reference, which can be accomplished 
by active stabilization of distributed interferometers 0, 
two-way roundtrip communication JjJ, or one-way auto- 



compensation using photon pairs [(| . The second recently 
proposed strategy exploits the advances in the theory of 
noiseless subsystems [t| to construct joint polarization 
states of several photons that are immune to collective 
depolarization (3, 111 liij • From the fundamental perspec- 
tive, such states enable quantum communication without 
a shared reference frame . Most of the work in this 
direction is based on results obtained for an ensemble of 
qubits that undergo an identical symmetrized interaction 
with the environment. In the optical domain, this model 
translates into a sequence of pulses containing one photon 
each and prepared in joint, usually entangled, polariza- 
tion states. 

Sequences of identically polarized pulses employed in 
phase encoding, or trains of single-photon pulses used 
to implement noiseless subsystems, form very particular 
classes of states in the entire Hilbert space describing a 
set of optical modes travelling in a fiber. The central aim 
of this paper is to analyze the structure of the complete 
Hilbert space under collective depolarization and to de- 
termine its potential for implementing quantum commu- 
nication protocols. This approach will provide a general 
treatment of depolarization in optical fibers that includes 
a number of previously proposed schemes as special cases. 
Furthermore, our analysis will reveal the existence of hy- 
brid noiseless subsystems that combine both polarization 
and multiphoton phase encoding in a non-trivial way. 
These hybrid subsystems are optimal in terms of required 
resources, namely the total numbers of temporal slots 
and photons. 

We shall consider general, multiphoton states of optical 
radiation prepared over N temporal slots that experience 
identical birefringence. The polarization transformation 
between a pair of annihilation operators describing two 
orthogonally polarized modes that occupy one temporal 
slot is given by an element ft of the Lie group U(2). 
The element J~i can be decomposed into a product of the 
overall phase factor, which we will denote by e _8Q ( n ) , and 
the remaining matrix e la ^£l with the determinant equal 
to one which belongs to SU(2). The element O induces a 
unitary transformation U(Sl) in the corresponding two- 
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mode Fock space. For our purposes it is natural to use 
the Schwinger representation |l2| in which the Hilbert 
space is decomposed into subspaces with a fixed number 
of photons, which we will denote by I. The polarization 
transformation in the /-photon subspace is given by (I + 
l)-dimensional representation of the rotation group for 
spin 1/2. Explicitly, the decomposition of U(il) in the 
Schwinger representation has the form: 



j/(n) =0e- iiQ < n >P'/ 2 (e iQ(n) «) 



(1) 



where the matrix V l ^ 2 (e la ^fl) is the corresponding ele- 
ment of the (7 + l)-dimensional irreducible representation 
of SU(2). 

The transformation due to depolarization of the entire 
quantum state g of radiation contained in TV temporal 
slots is given by the integral: 
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where d£l is the invariant Haar measure in U(2). We are 
now interested in identifying the subspaces of the entire 
Hilbert space that preserve quantum coherence. For this 
purpose, we need to expand the TV-fold tensor product 
of the unitary transformation U(fl) into a direct sum of 
irreducible representations. If we first divide the TV-slot 
Hilbert space into sectors containing in total exactly L 
photons, this expansion will in general take the following 
form: 



®N 



0< 

L=0 



-iLa(Sl) 



L/2 

j=(L mod 2)/2 

(3) 

In the above expression, the integers K J NL define how 
many times the spin-j representation appears in the sec- 
tor of L photons distributed between TV slots, and we 
have symbolically denoted by a dot the operation of tak- 
ing a direct sum of K 3 NL replicas of this representation. 

If we now insert the decomposition given in Eq. Q 
into Eq. J2Jl, the invariant integration over CI will com- 
pletely remove coherence between subspaces with differ- 
ent L, and also between any two sectors with different 
j. The latter fact is a result of the canonical property of 
the products of rotation matrix elements integrated over 
SU(2). Consequently, the sector of the Hilbert space with 
given L and j can be treated from the point of view of 
quantum information processing applications as a noise- 
less subsystem with dimensionality K J NL . This system 
can be used for example to implement quantum cryptog- 
raphy protocols based on qudits |l3l| , or any other quan- 
tum communication task. Another application for which 
the above decomposition is relevant is the transmission 
of classical information 0, ^| . Then the system com- 
prising TV slots and containing at most L photons can be 



used to encode £l'=o T,jJ(l> mod 2)/2 K nl> distinguish- 
able classical messages. 

Let us now find and discuss explicit values of multi- 
plicities K J NL . First, we will analyze a restricted case 
when at most only one photon is allowed to occupy a 
temporal slot. We can then use the well-known result 
that for an ensemble of L qubits the allowed values of 
j are (L mod 2)/2, . . . , L/2, and that the spin-j algebra 



appears 



2j+l ( L+l 
L+l \L/2-j) 



times in a direct-sum decomposi- 
tion 0| . When the qubits are realized as polarization 
states of single photons and we have additional freedom 
to distribute the photons between TV temporal slots, these 
values are multiplied by the factor (^) which gives the 
number of configurations to occupy L out of TV slots. 
Thus the dimensionality K J NL of a spin-j noiseless sub- 
system (labelled with a tilde to distinguish it from the 
general scenario with no constraints on the number of 
photons in a slot) is given by 
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TV 
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L) L + l \L/2-j 



(4) 



The dimensionalities of noiseless subsystems for TV, L < 6 
are collected in Tabled It is instructive to compare these 
numbers with the case of pure phase encoding, when all 
the input photons are prepared in an identical polar- 
ization. This implies that their joint state of polariza- 
tion belongs to the completely symmetric subspace with 
the highest value of j, equal to L/2. Therefore phase- 
encoding schemes are included in Table [I] as entries with 
j = L/2, and for clarity they have been underlined. An 
interesting question is whether the dimension of these 
subspaces can be enhanced by preparing input photons 
in non-trivial polarization states. The answer is positive 
for TV > 4 and given by the highest entry in each col- 
umn marked with an asterisk. These entries correspond 
to hybrid encodings, where both the phase and the po- 
larization are exploited to protect quantum coherence. 
The first non-trivial hybrid noiseless subsystem occurs 
for three photons in four temporal slots, and it is ob- 
tained by extending the j = 1/2 subspace of three qubits 
by an additional empty temporal slot. In general, a sim- 
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pie analysis of the sign of the difference K 
shows that for given TV and L the largest subsystem is 
obtained for the highest integer (L even) or half-integer 
(L odd) value of j satisfying j < \JL + 2/2. Therefore, 
the largest noiseless subsystem will usually have a hybrid 
character. 

It is instructive to identify some earlier proposals in 
the literature with entries in Table [I] For example, the 
two protocols for quantum key distribution proposed in 
Ref. |s| are based on the cases TV = L = 4 and TV = L = 3 
with the corresponding values of j equal respectively to 
and 1/2, and two-dimensional noiseless subsystems. The 
essential advantage of these protocols is that they do 
not require a phase reference, in contrast to the most 
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TABLE I: Multiplicities K 3 NL of spin-j representations for L 
photons distributed between N slots with at most one pho- 
ton occupying each slot. Asterisked entries denote hybrid 
noiseless subsystems with the highest dimensionality in each 
column. 
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straightforward phase encoding scheme N = 2, L = 1 
that gives a noiseless subsystem of the same dimension. 
Similarly, the alignment-free test of Bell's inequalities 
proposed by Cabello @ uses N = L = 4, and Table U 
suggests that a similar scheme should be possible also 
for N = L = 3. 

Given present advances in the area of scalable linear- 
optics quantum computing it is interesting to investigate 
the asymptotics of hybrid noiseless subsystems for large 
N. Two relevant figures of merit are the average num- 
ber of logical qubits per slot that can be encoded in a 
sequence of N temporal slots, given by the expression: 

4 Q) = ilog 2 ^max^ (5) 

and the classical communication capacity of a train of N 
slots renormalized by the number of slots: 

4 c) = ^io S2 £ £ &nl • (6) 

\L=0j=(L mod 2)/2 J 

We found numerically by calculating these two quanti- 
ties up to N — 10 4 that they tend to the same con- 
stant value log 2 3, with asymptotics following a power 
law with the approximate value of the exponent C^' C ^ ~ 
log 2 3— 0(iV~ - 84 ), as shown in Fig.GJ In the same graph, 
we also depict the average number of photons per slot 
that are required to achieve these capacities, with the 
same asymptotic limit of (L) /N = 2/3. Thus we see that 




N 

FIG. 1: Asymptotics of the capacities C ® (•) and (■) 

with the corresponding average numbers of photons {L)/N 
for the case of quantum (o) and classical (□) communication. 

asymptotically the system behaves as if composed of a 
train of N noiseless qutrits, each spanned by the vacuum 
state and two one-photon states with orthogonal polar- 
izations. 

Let us now relax the constraint of having at most one 
photon per temporal slot and analyze the general sce- 
nario in which photons can be arbitrarily allocated in 
the temporal slots. We will derive a closed recursion for- 
mula for the multiplicities K J NL . First, let us simplify 
the notation by rewriting Eq. J2J symbolically in terms 
of the corresponding algebras of the noise operators: 

L/2 

A NL = K* NL -At (7) 

j=(L mod 2)/2 

Here Anl is the algebra for L photons in N slots, and 
the Ai denote standard irreducible spin-j algebras. The 
recursion formula will be based on a relation linking 
the representation multiplicities K 3 NL for L photons dis- 
tributed between N slots to those describing the system 
with the number of slots reduced by one. When we dis- 
tribute L photons between N slots, we can put an arbi- 
trary number L' < L of them into N — 1 slots we had 
before, and the remaining L — L' photons into the ad- 
ditional slot. Obviously, the polarization transformation 
of these L — L' photons in the new slot is governed by 
the algebra A^ L ~ L Assuming that the decomposition 
of the algebra An-i,l' f° r the photons in the remaining 
N — 1 slots is known in the form of Eq. J7J , this gives us 
the following relation: 

a nl =®awJ /v; a). 

L'=0 \j>=(L' mod 2)/2 / 

(8) 

We can now use the standard decomposition of the tensor 
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product of spin algebras: 

(L-L')/2+f 

A {L - L 'V 2 ® A = A j (9) 

j=\(L-L>)/2-j>\ 

to obtain an expression for Anl in the form of a triple 
direct sum: 

L L'/2 (L-L')/2+f 

Anl=® K^A* 

L'=Oj'=(L' mod 2)/2j=\(L-L')/2-j'\ 

(10) 

We now need to rearrange the summations in order to 
bring the above expression to the form of Eq. (JJJ . This 
can be done by replacing the first two summation indices 
L' and j' by a pair of new integer parameters /i = L' /2+ j' 
and v = L' /2—j'. Then the two-dimensional grid (L 1 , j') 
of the summation points with L' — 0, 1, . . . , L and j' = 
(£' mod 2)/2, (£' mod 2)/2 + 1,.. . , L'/2 can be written 
as {[i + v, ^-^-) with the new parameters running through 
ranges v = 0, 1, . . . , /2J and /U = z/, ^ + 1, . . . , L— v. In 
the new parametrization, Eq. (|10(l takes the form: 

L-L/2J t-u L/2-u 

*n L = e © ^ri^^- (ii) 

If we now want to perform first the summation over j, 
we need to identify the allowed ranges of /i and v as a 
function of j. Considering the limits of the sum over j in 
Eq. (|llfl . we obtain constraints on /x and ^ in the form 
v < L/2 - j and L/2 - j < fi < L/2 + j. These con- 
straints are always stronger than the original limits of the 
sums over \x and v in Eq. 111(1 . and they can replace the 
latter. Consequently, swapping the order of summations 
gives the following recursion formula for the multiplicities 

L/2-j L/2+j 
K j NL= E E K{ N-^l- (12) 

In Table HT1 we calculate explicitly the values of -ft^x for 
TV < 8 and L < 4. As before, for each pair of ./V and L the 
underlined entry with the highest j = L/2 corresponds to 
pure phase encoding which can be used as a benchmark 
to indicate whether exploiting the polarization degree of 
freedom can yield a larger noiseless subsystem. The opti- 
mal noiseless subsystems for given N and L in most cases 
turn out to be hybrid, as shown by asterisked entries in 
Table CH 

In conclusion, we have analyzed the structure of the 
Hilbert space for optical radiation contained in an arbi- 
trary number of temporal modes and undergoing a collec- 
tive depolarization, which is one of the dominant dccoher- 
ence mechanisms in optical fibers. We have demonstrated 
that optimizing communication capacity requires in gen- 
eral a combination of phase and polarization encodings 



TABLE II: Multiplicities K 3 NL for L photons distributed over 
slots, including arbitrary multiphoton states in each tem- 
poral slot. Underlined entries correspond to pure phase- 
encoding schemes, whereas asterisks mark optimal hybrid 
noiseless subsystems for given N and L. 
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into joint hybrid schemes. Although practical implemen- 
tation of these schemes may require complex preparation 
and measurement procedures, present experimental effort 
towards the realization of scalable linear-optics quantum 
computing allows for optimism that the necessary tech- 
niques will become available. As the hybrid noiseless 
subsystems utilize relative phases between pulses in a 
temporal sequence, the communicating parties require in 
general a shared phase reference. An interesting question 
is whether for hybrid subsystems this requirement can be 
circumvented in a way similar to the autocompensation 
technique introduced by Walton et al. 
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